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Observation of the GZK Cutoff Using the HiRes Detector 

D. R. Bergman a * presented on behalf of the High Resolution Fly's Eye Collaboration 

a Rutgers, The State University of New Jersey 
Department of Physics and Astronomy 
Piscataway, New Jersey, USA 08854 

The High Resolution Fly's Eye (HiRes) experiment has observed the GZK cutoff. HiRes observes two features 
in the ultra-high energy cosmic ray (UHECR) flux spectrum: the Ankle at an energy of 4 x 10 18 eV and a high 
energy suppression at 6 x 10 19 eV. The later feature is at exactly the right energy for the GZK cutoff according to 
the E-l/2 criterion. HiRes cannot claim to observe a third feature at lower energies, the Second Knee. The HiRes 
monocular spectra are presented, along with data demonstrating our control and understanding of systematic 
uncertainties affecting the energy and flux measurements. 



HiRes has observed the GZK Cutoff. The GZK 
cutoff was predicted 40 years ago, Greisen Q and 
Zatsepin and Kuzmin |2| predicted the existence 
of a sharp reduction in the flux of ultra-high en- 
ergy cosmic rays (UHECR's) at energies above 
about 6 x 10 19 due to photopion production inter- 
actions between UHECR protons and the cosmic 
microwave background radiation (CMBR). 

Since the prediction of the GZK cutoff was 
made, a number of experiments claim to have 
observed events with measured energies above 
1 x 10 20 eV. Volcano Ranch 3j observed one event 
in 1963, even before the GZK Cutoff was pro- 
posed. However, the knowledge of the expected 
lateral distribution function to use in estimating 
the total number of shower particle was lacking 
at that time. SUGAR claimed three events 4] 
and Haverah Park observed four0 . SUGAR used 
only muon detectors, however, and had problems 
with after pulsing in its PMT's. One should "be 
cautious about taking the energies ascribed to 
the Sydney events ... as providing definitive ev- 
idence against a cutoff in the cosmic ray energy 
spectrum0 i '' The Haverah Park data was later 
reanalyzed |2j using Corsika to establish the rela- 
tion between /?600 and the energy; all the high 
energy events were reanalyzed to lower energies, 
below the GZK cutoff. More recently, both the 
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Yakutsk ArrayQ and AGASA0] have claimed to 
events above 10 20 eV, but Yakutsk claims that 
it's events are consistent with what is expected 
from the GZK process while AGASA had claimed 
otherwise, but is now reconsidering their energy 
scale 

The point of going through this short history 
is to note that it is hard to measure the energy 
well with a surface array (which all of the above 
are), because one is never looking at the bulk of 
the shower. This makes the energy determina- 
tion subject to the poorly understood systematics 
of shower modeling, and it is necessary to mea- 
sure the energy well to observe a break in such 
a steeply falling spectrum, ft is easier, on the 
other hand, to measure the energy well with a 
fluorescence detector, such as HiRes. Because a 
fluorescence detector does observe the bulk of the 
shower, one is left only with systematic uncertain- 
ties which are easier to understand and control. 

Two such systematic uncertainties are number 
of photons emitted by the shower as it passes 
through the atmosphere, the fluorescence yield, 
and the clarity of the atmosphere through which 
one observes those emitted photons. The fluo- 
rescence yield has been well studied in the lab, 
and the differences between three recent measure- 
ments is only at the 6% level [HE El. The at- 
mospheric uncertainty, on the other hand, can be 
controlled by going to a site with a clean stable 
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atmosphere, like Dugway, Utah, and then by mea- 
suring the scattering properties of the air with 
lasers as we have done in HiRes[l^|. 

A third systematic uncertainty for a fluores- 
cence detector is the size of the aperture as a 
function of energy. Traditionally, ground arrays 
were thought to have constant aperture above the 
energy where the array becomes 100% efficient. 
In fact, given energy and core location resolu- 
tions this is not quite true. Fluorescence detec- 
tors, however, have a manifestly changing aper- 
ture because higher energy showers are brighter 
and can be seen farther away. The problem then 
becomes one of determining the size of the aper- 
ture at a given energy by a Monte Carlo simula- 
tion. However, this is not a case where we have 
replaced one poorly understood simulation, ex- 
tensive air shower development, which is the bane 
of ground array experiments, with another poorly 
understood simulation. The aperture calculation 
for a fluorescence detector uses only well under- 
stood physics and, more importantly, is amenable 
to detailed verification through the technique of 
Data/MC Comparisons. It is just these detailed 
comparisons that comprise the bulk of this talk 
and give us the great confidence we have in our 
aperture calculation. 

I will begin this discussion of the aperture cal- 
culation by noting that the aperture is, in ac- 
celerator physics parlance, the acceptance of the 
detector, and that acceptances are routinely cal- 
culated by MC simulation: 

= Acpt(E')Area T r2 T (2) 

where E' is the reconstructed energy, while E is 
the thrown energy. Nt and Nr are the num- 
bers of events generated ("thrown") and recon- 
structed, respectively. For the calculation to be 
accurate, one must put the right distribution of 
showers in the right places (to get the area) , at the 
right angles (to the right solid angle), and with 
the right amount of light (to get the right energy 
scale). What is more, one must put in the right 
fluctuations in signal for the difference between 
reconstructed energy (E') and thrown energy (E) 
to be significant. Fortunately, these distributions 



are easy to check in Data/MC Comparisons. 

The first comparisons are for the distances 
of the showers from both HiRes-I and HiRes-II 
(events in the right places). These are shown in 
Figure d for the distance to the impact point in 
HiRes-I in three energy bins, and in Figure [3 for 
the impact parameter for showers from HiRes- 
II. The second set of comparisons shows the 
distributions of measured angles (events at the 
right angles): the zenith angle for HiRcs-I events 
in Figure 02 the angle between the shower and 
the ground in the shower-detector plane, tp, for 
HiRes-II events in Figure 01 A third comparison, 
in Figure El shows the brightness distribution for 
events at HiRes-II (events with the right amount 
of light). The final comparison, in Figure|Sl shows 
the distribution of x 2 for a fit of the times of tubes 
at HiRes-II to a linear function of time vs. angle 
(for resolution). Since the time .vs angle relation 
should be curved, the x 2 distribution has a large 
tail. 

Having built confidence in our aperture calcu- 
lation, we move on to the data measured by HiRcs 
in monocular mode. This is shown in Figured It 
is this distribution which we will be fitting later 
to determine the significance of our observation of 
the GZK Cutoff. The calculated exposure (aper- 
ture times time) is shown in Figure |H1 Dividing 
the event distribution by the exposure gives the 
spectrum, which is shown in Figure El 

We performed a number of systematic cross 
checks on these data to ensure their accuracy. 
First wc limited the distance to showers to 10 and 
15 km. We saw no systematic change in the spec- 
tra in these cases, only variations due to reduced 
statistics. What is more, the Ankle feature is vis- 
ible in the event distribution in the 10 km case. 
Second, we made a spectrum from the events 
which are in both the HiRes-I and HiRes-II sam- 
ples, using the precision shower geometry mea- 
surement available in this case. All the other de- 
tails of the analysis are the same as in the monoc- 
ular analyses, and the spectrum agrees with our 
monocular spectra. This shows that the energy 
resolution in monocular mode is good enough not 
to distort the spectrum from the stereo measure- 
ment beyond the statistical power of the data we 
have collected. This was not true of Fly's Eye. 
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Figure 1. A Data/MC Comparison of the dis- 
tance to the impact point from HiRes-I. The three 
panels show events in the different energy ranges 
as labeled. Data is shown as points with error 
bars, the MC as a histogram. 



Figure 2. A Data/MC Comparison of the impact 
parameter, Rp, from HiRcsTI. The top panel 
shows the data as points with error bars, the MC 
as a histogram. The MC has about five times the 
statistics as the data. The bottom panel shows 
the ratio of data to MC, with a linear fit over bins 
with more than 10 data events. 



We also tested the aperture at 35 km, by recon- 
structing vertical laser shots at this distance. A 
variety of laser intensities were used, which we 
linked to cosmic ray energies based on the bright- 
ness of the tracks. Our detector becomes fully 
efficient for vertical showers at 35 km for showers 
with energies above 60 EeV. 

Because we use a fit the Gaisser-Hillas profile 
to determine the shower energy, we compared the 
average shower profile in data to that calculated 
in MC using the Gaisser-Hillas profile. In each 
case the showers are normalized by A max (using 
the age parameter) and A max . Then the actual 
shower measurements are averaged in bins of the 
age parameter. The comparison is shown in Fig- 
ure This shows that the Gaisser-Hillas profile 
is appropriate to use for measuring shower ener- 
gies. 

To check the systematic uncertainty from 



shower modeling, we calculated the aperture us- 
ing both QGSJet and Sibyll[l5j. It is important 
to note that we adjust the proton-to-iron ratio 
in each case to give identical mean Xmax values 
as a function of energy. With this requirement, 
we calculate identical apertures for both models 
within the MC statistics (5x the data) available. 

We cannot claim to observe the Second Knee 
in the HiRes-II spectrum. The data statistics 
are limited and the aperture has a larger un- 
certainty than at energies above 10 18 eV due to 
composition uncertainties and their effect on the 
aperture^]. Because of this, we will only be fit- 
ting the spectra above 10 17 5 eV. 

We now move on to broken power law fits of the 
two HiRes monocular spectra. We will be fitting 
the spectra using a binned maximum likelihood 
method which uses the number of observed events 
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Figure 3. A Data/MC Comparison of the zenith 
angle in HiRes-I. The three panels show events 
in the different energy ranges as labeled. Data 
is shown as points with error bars, the MC as a 
histogram. 



Figure 4. A Data/MC Comparison of the shower- 
detector plane angle, ip, from HiRes-II. The top 
panel shows the data as points with error bars, 
the MC as a histogram. The MC has about five 
times the statistics as the data. The bottom panel 
shows the ratio of data to MC, with a linear fit 
over bins with more than 10 data events. 



in a bin and the number of events expected. This 
will allow us to use empty bins in the data where 
there is a significant number expected from the 
power law. It also allows us use both data sam- 
ples independently A correction will be made at 
the end for events seen by both detectors. The 
number of expected events is given by the broken 
power law flux times the measured exposure. We 
will begin with no break points and add floating 
break points until there is no improvement in the 

x 2 - 

A fit of both datasets to a simple power law 
gives a \ 2 of 162 for 39 degrees of freedom, with 
a spectral index 7of3.13±0.01. This is clearly 
not a good fit. Adding one floating break point 
gives a much better fit, and the break point finds 
the position of the Ankle. The \ 2 1S now 68 for 
37 degrees of freedom, with the break point at 
4.3 ± 0.5 EeV, a spectral slope below the Ankle 



of 7i = 3.24 ± 0.02 and a spectral slope above 
the Ankle of 72 = 2.89 ± 0.03. Adding a third 
break point reduces the \ 2 by 33 to 34.7 for 35 
degrees of freedom. The breakpoints are at 4.5 
EeV and 56 EeV, while the three spectral indices 
are 3.24 ± 0.02, 2.81 ± 0.03 and 5.4 ± 0.7. This fit 
is shown in Figure Adding a third breakpoint 
does not significantly reduce the \ 2 ■ 

We calculate the significance of the break at 56 
EeV by comparing the number of events expected 
above 10 19 8 eV if the spectral index continued at 
2.81 to the actual number of events observed. In 
this case we expect 51.1 events but observe only 
15. The Poisson probability of observing 15 or 
fewer events when expecting 51.1 is 2.9 x 10 -9 . 
In actuality, we have counted one event observed 
in stereo twice in our 15 observed events and have 
also counted expected events from overlapping ex- 
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Figure 5. A Data/MC Comparison of the event 
brightness, total signal divided by track length, 
for HiRcs-II. The top panel shows the data as 
points with error bars, the MC as a histogram. 
The MC has about five times the statistics as the 
data. The bottom panel shows the ratio of data 
to MC, with a linear fit over bins with more than 
10 data events. 



Figure 6. A Data/MC Comparison of the dis- 
tribution of x 2 /rmDOF for a linear fit of tube 
time vs angle, for HiResTI. The top panel shows 
the data as points with error bars, the MC as 
a histogram. The MC has about five times the 
statistics as the data. The bottom panel shows 
the ratio of data to MC, with a linear fit over bins 
with more than 10 data events. 



posure from the two sites. When each of these is 
taken into account we get a revised expectation of 
44.9 events while observing 14. The Poisson prob- 
ability in this case is 7 x 10~ 8 . For comparison 
the fraction of the area in one tail of a Gaussian 
beyond 5cr is 3 x 10~ 7 and beyond 6<7 is 1 x 10~ 9 . 
Thus we claim a significance of over 5ct for our 
measurement of a high energy break point in the 
spectrum. 

One way to compare the energy of this break- 
point to that expected of the GZK Cutoff is to 
use the E 1 / 2 method of Berezinskv[l^|. E 1 / 2 is 
the energy at which the integral spectrum falls to 
half of what would be expected in the absence of 
the feature. This HiRes monocular spectra are 
displayed as integral spectra in Figure ^] along 
with the two break point fit mentioned above and 



the extension of this fit used to measure the sig- 
nificance of the high energy break point. Taking 
the ratio of the measured points to the extended 
fit allows us to measure Ei/ 2 - This is shown in 
FigurelEl We find log£ 1/2 = 19.73 ± 0.07. The 
prediction Berezinsky et aZ|l6| for the GZK effect 
for a wide range of spectral indices is log Exp = 
19.72. In this way we have linked the energy of 
our measurement of a high energy break point to 
the GZK Cutoff. 

In conclusion, HiRes has made a measurement 
of the UHECR spectrum using its two detectors 
in monocular mode. The aperture used for this 
measurement is well understood as verified by 
Data/MC comparisons. We observe a high en- 
ergy break in the spectrum at the energy to be 
the GZK Cutoff and with a significance of over 
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Figure 7. The number of events collected in each 
energy bin in monocular mode by HiResT and by 
HiRcs-II. The jump in HiRes-II statistics at 10 19 
eV is due to the bins being twice as big in logE. 



Figure 8. The calculated exposure of HiRes-1 and 
HiRcs-II in monocular mode. The exposure has 
been divided by E 2 to make direct comparisons 
to the event distribution easier. 



5(7. Thus we claim an observation of the GZK 
Cutoff. 

This work is supported by US NSF grants 
PHY-9321949, PHY-9322298, PHY-9904048, 
PHY-9974537, PHY-0098826, PHY-0140688, 
PHY-0245428, PHY-0305516, PHY-030098, and 
by the DOE grant FG03-92ER40732.We grate- 
fully acknowledge the contributions from the 
technical staffs of our home institutions. The co- 
operation of Colonels E. Fischer and G. Harter, 
the US Army, and the Dugway Proving Ground 
staff is greatly appreciated. 

References 

[1] K. Greisen, Phys. Rev. Lett. 16 (1966) 748- 
750. 

[2] G. T. Zatsepin, V. A. Kuz'min, J. Exp. 
Theor. Phys. Lett. 4 (1966) 78-80, [ZhETF 
Pis'ma 4 (1966) 114-117]. 



[3] J. Linsley, Phys. Rev. Lett. 10 (1963) 146- 
148. 

[4] M. M. Winn, et al., J. Phys. G 12 (1986) 
653-674. 

[5] M. A. Lawrence, R. J. O. Reid, A. A. Wat- 
son, J. Phys. G 17 (1991) 733-757. 

[6] M. Nagano, A. A. Watson, Rev. Mod. Phys. 
72 (2000) 689-732. 

[7] M. Ave, et al., Astropart. Phys. 19 (2003) 
47-60, [astro-ph/0112253]. 

[8] V. P. Egorova, et al., in Proceedings of the 
Cosmic Ray Internation Seminar, Catania 
2004, Nucl. Phys. B Proc. Supl. 136 (2004) 
3-11. 

[9] M. Takeda, et al., in Proceedings of the 28th 
ICRC (Tsukuba), pp. 381-384, 2003. 

[10] M. Teshima, at this conference. 



7 




c 

1 

0.8 
0.6 
0.4 
0.2 




DA1A-Vor.teCa<>o 















* 


>■ 

I 








• 


I 










t 


> 

• 








1 DATA 


t 

■ 






1 




i 


W 



0/. 0.6 0,6 



1.4 1.6 



20 20.5 

log.„(E) (eV) 



Figure 10. A Data/MC comparison of the average 
profile for showers in HiRes-II. 



Figure 9. The measured UHECR spectra from 
HiResT and HiRes-II in monocular mode. Also 
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son. The HiResT and HiRes-II spectra have been 
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Figure 11. The integral spectra measured by 
HiResT and HiRes-II. Also shown are the inte- 
grated broken power law fits to the differential 
spectra. The black line shows the fit with two 
break points; the red line shows the same fit in 
which the high energy break point has been re- 
moved with no change in the other parameters. 
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Figure 12. The ratio of the HiRes-I and HiRes-II 
integral spectra to the expected integral flux in 
the absence of a GZK cutoff. Only the HiRes-1 
values are used to make an estimate of E 1 / 2 by 
interpolating between the central value and lcr 
uncertainty limits. 



